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INTRODUCTION 

The Numerical Aerodynamic Simulator (NAS) 
epitomizes the coming of age of supercomputing and 
opens exciting new horizons in the world of numerical 
simulation. The technological benefits likely to result 
from the use of the NAS may surprise even its most 
optimistic advocates. All this makes me both pleased 
and honored to be a participant of this conference 
which celebrates the inauguration of the NAS. 

In this article I will give an overview of supercom- 
puting at Lockheed Corporation in the area of Compu- 
tational Fluid Dynamics (CFD). This overview will fo- 
cus on developments and applications of CFD as an 
aircraft design tool and will attempt to present an as- 
sessment, within this context, of the state-of-the-art in 
CFD methodology. Although CFD is being actively de- 
veloped and applied throughout Lockheed Corpora- 
tion, the material of this paper draws heavily from ac- 
tivities that have taken place or are underway at 
Lockheed-California Company. 

Of course, supercomputing in the aerospace indus- 
try is not limited to CFD. Many other disciplines either 
are already, or have the potential of, benefiting from 
supercomputing: structural analysis and optimization, 
vehicle signature (electromagnetic, acoustic, thermal) 
prediction, and signal processing, to name the obvious 
ones. In other areas, where the supercomputing revolu- 
tion has not yet made its full impact, the payoffs are 
likely to be extraordinary. One such area is manufac- 
turing. Supercomputers will make possible the effective 
integration of all phases and aspects of complex manu- 
facturing processes and are bound to become essential 
in the development and operation of a new generation 
of sophisticated, versatile, and 'intelligent' robotic sys- 
tems. 

In spite of the wide applicability of supercomput- 
ing, I have chosen to limit my discussion to CFD and 
its application to aircraft design because of the follow- 
ing factors: 1) time and space constraints; 2) personal 
bias due to my familiarity with the subject; and 3) CFD 
was one of the first disciplines to push the limits of 
'ordinary' or 'general purpose' computers, thus giving 
rise to the need for supercomputer capability. More- 
over, the practical feasibility of CFD application de- 
mands supercomputing power, and most future devel- 
opments in this area hinge on the availability of even 
greater computer speed and memory capability. 


SUPERCOMPUTERS 
AT LOCKHEED CORPORATION 

The term supercomputer is a relative one; what we 
consider supercomputer performance today may be- 
come 'ordinary' or 'general purpose' computer capa- 
bility in a few years. For the purpose of the present 
discussion, the term 'supercomputer' will denote a ma- 
chine with a performance capability of at least class VI 
in accordance with the following classification defined 
by the U.S. Department of Energy: 

Class III Sustained operating speeds between 0.6 and 
2 million floating point operations per sec- 
ond (megaflops) 

Class IV 2 to 6 megaflops 

Class V 6 to 20 megaflops 

Class VI 20 to 60 megaflops 

Class VII Greater than 60 megaflops 

Presently two supercomputers are in operation at 
Lockheed Corporation. The first one, a Cray IS, was 
installed at the Sunnyvale, California, plant of Lock- 
heed Missiles and Space Company (LMSC) in 1982. 
This supercomputer is basically dedicated to serving 
the needs of the various LMSC divisions. The second 
one, a Cray X-MP/24, became operational at the Kelly 
Johnson Research Center, Rye Canyon, California, in 
April 1986. 

The Cray X-MP/24 is a dual processor machine 
with an in-core memory capacity of 4 million 64 -bit 
words. It is front-ended by two Digital Equipment Cor- 
poration (DEC) computers in the Lockheed installation: 
a DEC 8600 for unclassified processing, and a DEC 
8300 for classified mode operation. This supercompu- 
ter, although it can be accessed remotely from all major 
Lockheed divisions, is primarily intended to cater to the 
Lockheed Aeronautical Systems Group companies. 

The principal motivation for the acquisition of 
supercomputer hardware was the perceived need for 
processing speed. The above supercomputers have 
been able to fulfill that need fairly well for most present 
applications. At the same time, they have offered the 
possibility of more advanced computations. As we gain 
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more experience with these computations, and as we 
attempt to tackle more demanding problems, the need 
for more computer memory becomes acutely evident. 
This is particularly true in the engineering applications 
of CFD. As will be indicated later on, the use of ad- 
vanced CFD methods in aircraft design applications 
requires fast-access memory capacity ranging from 16 
megawords upwards. 

THE CHALLENGE OF 
DESIGN APPLICATIONS 

The paper is mainly concerned with CFD develop- 
ments and applications within the context of aircraft 
design. Therefore, although it may appear obvious to 
many, it is important to emphasize the principal re- 
quirements that a CFD code ought to meet if it is to be 
used as an aircraft design tool, namely: 

1) It must be a predictive tool, i.e., the accuracy 
and reliability of the computation should not 
depend on too much fine tuning or tailoring of 
input parameters. In many advanced applica- 
tions, applicable experimental data are not 
available to "'calibrate * the results of numerical 
computations. 

2) It must be able to handle complex geometric 
configurations, including three-dimensionality, 
intricate aircraft components, and complete air- 
craft arrangements. 

3) It must yield useful results in a reasonably short 
period of time, consistent with usually pressing 
program schedules. 

4) Its application, although it may require a rea- 
sonable degree of expertise and experience, 
must not demand specialization to a degree that 
only its developers can effectively use it. 

The above requirements are extremely important if 
CFD is to realize its full potential as an effective design 
tool. As has already been pointed out in Reference 1, 
these requirements directly relate to the equation: 

EFFECTIVENESS = QUALITY X ACCEPTANCE 

This expression has no actual quantitative mean- 
ing, it is merely a symbolism to emphasize the axiom 
that the effect a given process has on the activity for 
which it is intended depends not only on the goodness 
of the process itself, but also on how widely it is used 
or accepted. In the present context, the first couple of 
items listed above impact the quality factor; the last two 
the acceptance factor. 


The ability to meet the above requirements truly 
differentiate between a research CFD method and one 
that is an engineering and design tool: for design appli- 
cations, a CFD method must meet all above require- 
ments with a reasonable degree of completeness. It is 
not uncommon to confuse research feasibility with de- 
sign application readiness. It is the professional respon- 
sibility of the CFD code developer to fully acknowledge 
this difference and make it clear to the non-expert. 

Unfortunately, there is a tendency amongst CFD re- 
searchers to greatly underestimate, or completely disre- 
gard, the important effort of turning a research, or pilot, 
code into a practical engineering tool. 

Thus, the principal challenge posed by design ap- 
plications to CFD practitioners is to convert what is 
feasible in principle (i.e., feasible in a research laborato- 
ry environment) into practical engineering tools that 
meet all four requirements listed at the beginning of 
this section. 

Some of the principal CFD codes for aircraft design 
application in use at Lockheed are described in the 
following sections. These codes have been singled out 
for discussion herein on the grounds that they have 
been or are being developed specifically to be practical 
engineering tools. As in most engineering disciplines, 
there is a hierarchy in CFD methodology, namely, the 
various computational codes can be classified accord- 
ing to levels of increasing complexity and accuracy. 
This also implies different requirements in terms of 
computational resources. The codes to be discussed, all 
of them in practical use at Lockheed, span the CFD 
methodological hierarchy. 


LINEAR CFD METHODOLOGY: PANEL 
METHODS 

Methods based on linear inviscid flow assumptions 
are more fully developed than nonlinear methods, and 
they are less demanding of computer resources. The 
most common linear methods are the so called panel 
methods. The majority of these methods are limited in 
the flow physics they can model: they cannot deal with 
shock waves, transonic flows, and separated flows. On 
the other hand, they are capable of modeling all geo- 
metric features of the vehicle, including the most intri- 
cate details. At Lockheed, the advanced low-order pan- 
el method QUADPAN for subsonic and supersonic 
flows is the principal tool for linearized analysis, Figure 
1. It has been under development and in use for a 
number of years. Reference 2. 
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Figure 1. Advanced low-order quadrilateral element panel method 
QUADPAN. 

In a panel method, the airplane to be modeled is 
represented by a large number of small contiguous sur- 
face panels. Usually these panels are flat and of a quad- 
rilateral shape. Elementary solutions of the linearized 
flow equation, known as sources and doublets, are as- 
signed to each panel. Figure 2. The individual strengths 
of these elementary solutions are determined by the 
application of the appropriate boundary conditions, 
e.g., no flow across solid surfaces. The shape of the 
distribution across a panel of these elementary solu- 
tions determine whether the method is low order (e.g., 
constant or linear variation) or high order (e.g., qua- 
dratic variation). 


(SOURCE) 



Figure 2. Concept of panel methods 


Until recently it was a common belief that a high 
order distribution was required for better accuracy. The 
use of high order distribution imposes a sizeable com- 
putational cost penalty. Special numerical techniques, 
discussed in reference 2, can greatly improve the accu- 
racy and robustness of low order singularity distribu- 
tions. QUADPAN's advanced low order approach em- 
bodies these techniques which give it an accuracy 
similar to that of the high-order methods but at a frac- 
tion of the computational cost. This has been verified 
in evaluations conducted by NASA and U.S. Air Force 
researchers, References 3 and 4. 


Figure 3 illustrates the result of a typical QUAD- 
PAN computation which makes evident the code's 
ability to model complex geometry. It shows a generic 
advanced fighter configuration with the corresponding 
elemental quadrilateral panels used to model the flow 
about it. The computed pressure distribution, at an an- 
gle of attack of 10 degrees at a mach number of 0.90, 
is shown in a color-coded display. Correlation of pre- 
dicted force and moment coefficients with wind tunnel 
data are presented in Figure 4. 



Figure 3. Surface pressure distribution computed by QUADPAN 
on generic fighter configuration. 





CD 


O QUADPAN 
EXPERIMENT 


Figure 4. Comparison of QUADPAN predictions with wind tunnel 
data for generic fighter model. 
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In spite of the inherent limitations of linearized 
theory, QUADPAN has proven to be a valuable design 
tool. Reference 5. This is due to its great geometric 
capability and to its ease of use and computational 
efficiency. These factors continue to make linearized 
methodology the workhorse for conceptual design ap- 
plication. Although the routine use of a linear method 
such as QUADPAN does not require supercomputer 
capability, the effectiveness of the method is greatly 
enhanced by it due to resulting reductions in computa- 
tional cost and cycle time. 

NONLINEAR CFD METHODOLOGY: EULER 
CODES 

Nonlinear methods have been undergoing inten- 
sive development. They either solve the full potential 
or the transonic small perturbation equations for invis- 
cid irrotational flows, or the Euler equations for in viscid 
rotational flows, or the Reynolds-averaged Navier- 
Stokes equations for viscous laminar and turbulent 
flows. Their increased accuracy and realism is attained 
at the cost of substantially larger computational re- 
quirements. Furthermore, the entire space surrounding 
the vehicle has to be covered with a computational 
mesh or grid in order to solve the equations of conser- 
vation of mass, momentum, and energy throughout the 
flow field as is required by the nonlinearity of the equa- 
tions. The generation of this computational grid for 
arbitrarily complex configurations is a difficult task, 
and it is of crucial importance for the practical applica- 
tion of these methods. 

For flow solution algorithms, the development of 
time-stepping, finite volume techniques is being em- 
phasized at Lockheed-Califomia Company. These 
techniques are essentially the same whether the Euler 
or the Navier-Stokes equations are being solved. The 
principal difference is the appearance, when solving 
the Navier-Stokes equations, of shear stress flux terms 
generated by the viscosity of the fluid and the turbu- 
lence of the flow. These flux terms are zero for the Euler 
equations, which model inviscid flows. 

Euler methods are now under advanced develop- 
ment and are beginning to see wide application in the 
design environment. Navier-Stokes methods are in an 
earlier stage of development and have only seen pio- 
neering applications. The effective application of these 
methods in a practical design environment makes the 
use of supercomputers mandatory. 

Although a number of different codes at Lockheed 
solve the Euler equations, the most advanced in terms 
of geometric capabilities is the TEAM (for 
Three-dimensional Euler Aerodynamic Method) code. 
Figure 5, the development of which is being partially 


funded by the U.S. Air Force Flight Dynamics Labora- 
tory, References 6, 7, and 8. The TEAM code is a modu- 
lar computational system consisting of preprocessors, 
grid generation routines, flow solver, and postproces- 
sors. Its flow solver is based on a variant of the finite- 
volume time-stepping algorithm proposed by Antony 
Jameson, et al. Reference 9. 

y t Iff Q dV + f f F.n dS = 0 FINITE VOLUME: 
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SUB/TRAN/SUPERSONIC + STRONG SHOCKS + VORTEX FLOWS 

Figure 5. Three-dimensional Euler Aerodynamic Method (TEAM) 
code. 


The ability of the TEAM code to handle arbitrarily 
complex geometries is due to its zonal structure. This 
means that the computational domain can be divided 
into several zones with interfacing boundary surfaces 
as shown in Figure 6. The various zones can have dif- 
ferent grid densities and the interfacing boundary sur- 
faces can accommodate different types of boundary 
conditions (fluid-to-fluid interface, solid surface, tran- 
spiration, free stream). All types of grid topologies (C- 
H, C-O, O-H, O-O, H-H, and combinations thereof) 
can be dealt with, Figure 7. Complete airplane configu- 
rations including wing, fuselage, canards, horizontal 
and vertical tails, flow-through nacelles and inlets, etc. 
are within the realm of the modeling capability of the 
TEAM code. 


• FACILITATES ANALYSIS OF REALISTIC AIRCRAFT 

• INCREASES COMPUTATIONAL EFFICIENCY 

• MORE ACCURATE FLOW SIMULATION 
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Figure 7. Typical computational mesh topologies. 


The code is able to handle subsonic, transonic, and 
supersonic flows. It can be applied to both steady and 
unsteady flows, and it automatically captures shock 
waves and vortex flows induced by sharp edges. The 
ability to capture vortex flows is of great practical sig- 
nificance, particularly in the design of combat aircraft. 
Leading edge vortices, like those appearing above a 
delta wing at high angle of attack, induce strong non- 
linear effects which cannot be predicted by classical 
potential flow theory. Euler methods now offer the 
capability of predicting these nonlinear effects for com- 
plex configurations. 


Figure 8 shows the results of a TEAM code compu- 
tation of the flow about a delta wing at high angles of 
attack. These results are compared with wind tunnel 
data and predictions made by classical linear theory. 
The improvement in accuracy provided by the Euler 
computations over linear theory is quite evident. Fur- 
thermore, the computed velocity field clearly displays 
the vortex-like structure which characterizes the flow 
about swept leading edges at high angles of attack. 
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Figure 8. Comparison of computations about delta wing with 
experimental data. 


The solution of the Euler equations appear to do an 
adequate job of predicting the development of vortex 
flows induced by sharp or highly swept edges, at least 
up to the point of vortex burst or breakdown. The geo- 
metrical flexibility provided by the TEAM code allows 
the analysis of complex vortex flow interactions about 
complete aircraft configurations. Our experience with 
the TEAM code indicates that in order to obtain reason- 
ably accurate results, fairly dense computational grids 
are needed. Typically, as many as 400,000 cells may be 
required for a configuration consisting of a wing, fuse- 
lage, and horizontal and vertical tails in a symmetric 
flight condition. 

Some Euler computations show the appearance of 
unsteadiness in vortical flows at conditions which coin- 
cide with the experimental observation of the onset of 
vortex breakdown. Other qualitative characteristics in- 
dicative of the effects of vortex burst have also been 
observed in these computations. An example of this is 
illustrated in Figure 9, which shows that the predicted 
effect of a leading edge strake upon vortex flow devel- 
opment correlates reasonably well with the wind tun- 
nel data. But more study and correlations are needed 
before ascertaining the usefulness of Euler solutions for 
the prediction of vortex burst. 


STRAKE OFF 
COMP. Cl = 0.82 
EXP. C L =0.89 



Afl = 0.5, a = 18° 


aCl: COMP= 0.32 
EXP = 0.34 

X/C «= 0.77 


STRAKE ON 
COMP. C L = 1.14 
EXP. C L = 1.23 






Figure 9. Effect of strake on leading edge vortex flow. 


The TEAM code is also being successfully used in 
high supersonic mach number applications. Figure 10 
presents the surface pressure distribution on one of the 
waverider configurations of Reference 10 at a mach 
number of 6 and at 4 degree angle of attack. Correla- 
tions of the computed pressure distributions with the 
available wind tunnel data are shown in Figure 1 1 . 

The value of Euler computations in design applica- 
tions have been amply demonstrated in many cases. To 
what extent they will be able to satisfy the most press- 
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Figure 10. Surface pressure distribution computed by TEAM 
code on waverider configuration at Mach number = 
6 and 4-degree angle of attack. 



LEGEND 
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........ AL = -4, TEAM © AL = -4.EXP 

AL = 4, TEAM □ AL = 4, EXP 

Figure 11. Comparison of predicted lower surface pressure dis- 
tribution with wind tunnel data. 


ing prediction requirements posed by the design of ad- 
vanced aircraft is far from being fully established. 
Much work, particularly in terms of advanced compu- 
tations and correlations with appropriate high quality 
experimental data, has to be done before the realm of 
practical applicability and validity of Euler solutions 
can be determined. 


NONLINEAR CFD METHODOLOGY: 
NAVIER-STOKES CODES 

The solution of the Navier-Stokes equations offer 
the potential of modeling all the physics of a viscous 
compressible fluid. The calculation of the viscous ef- 
fects from first principles for all levels of turbulence is 
presently beyond the realm of practical feasibility due 
to the lack of adequate computing power. Therefore, 
the Reynolds-averaged form of the equations is used in 
engineering applications. The Reynolds-averaged 
Navier-Stokes equations include the viscous terms but 
semi-empirical models must be introduced to represent 
the flow turbulence. 

A wide variety of Reynolds - averaged 
Navier-Stokes codes are in use at Lockheed, including 
codes with real gas effects and chemistry. But the most 
advanced from the point of view of geometric capabili- 
ty is the TRANSAM (for Three-dimensional Reynolds- 
Averaged Navier-Stokes Aerodynamic Method) code, 
Figure 12. Basically, it constitutes an extension of the 
modularized TEAM computational system to which 
momentum fluxes due to both viscous and Reynolds 
stresses have been added. Therefore, it possesses geo- 
metric flexibility and generality similar to those of the 
TEAM code, and it is equally able to analyze subsonic, 
transonic, and supersonic flows. 


4///QdV + // F.ndS= A // F v .n dS FINITE VOLUME: 
dt R N 



SUB/TRAN/SUPERSONIC + STRONG SHOCKS + VISCOSITY 


Figure 12. Three-dimensional Reynolds-averaged Navier-Stokes 
Aerodynamic Method (TRANSAM) code. 

In the present version of TRANSAM, either the full 
Reynolds-averaged or the thin shear layer approxima- 
tions to the Navier-Stokes equations can be solved at 
the user's option. The presently implemented turbu- 
lence model is an algebraic one, a modified version of 
the Baldwin-Lomax eddy viscosity model. Since it is 
architectured in a multiple zonal structure, it is possible 
to solve different equation types in each zone. For in- 
stance, the thin shear-layer approximation to the Reyn- 
olds-averaged Navier-Stokes equations can be solved 
in zones close to solid boundaries where boundary- 
layer behavior is to be expected; all the shear stress 
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terms can be accounted for in zones where fully sepa- 
rated flow is likely to occur; and, finally, the Euler 
equations can be used to model the flow for the remain- 
ing zones. This approach yields obvious economies in 
both computer execution time and memory require- 
ments. 

An example of the potential improvement in accu- 
racy that can be obtained with a Navier-Stokes compu- 
tation is given by the analysis of the supercritical RAE 
2822 airfoil shown in Figure 13. Both viscous and invis- 
cid calculations done with TRANSAM are compared 
with experimental data. The improvement that the vis- 
cous computation yields in the prediction of shock loca- 
tion and overall character of the pressure distribution 
is quite evident. The computed velocity profiles in the 
boundary layer, shown before and after the shock in 
Figure 14, show the correct behavior. 



Figure 13. Comparison of computed surface pressure distributions 
with wind tunnel data for RAE 2822 supercritical air- 
foil. 


We have not yet attempted the computation of sep- 
arated flows about aircraft at high angles of attack be- 
cause of lack of adequate memory capacity. Even using 
the zonal approach provided by TRANSAM, a mini- 
mum of 120 mega words of memory would be needed 
for such a computation. Experience with separated 
flows about simpler configurations indicates that the 
presently available algebraic turbulence models are in- 
adequate to predict strongly separated flows with rea- 
sonable accuracy and consistency. Whether more so- 
phisticated models (such as the two-equation or Reyn- 
olds stress models) will significantly improve the accu- 
racy of separated flow computations remains to be de- 
termined. 



U/UINF 


Figure 14. Before and after shock boundary layer profiles com- 
puted by TRANSAM code. 

In general, the various Navier-Stokes computa- 
tions of three-dimensional configurations that have 
been performed to date suffer from lack of adequate 
computational mesh resolution, and from the limita- 
tions of the turbulence models being used. Although 
some calculations that have recently been presented 
elsewhere look impressive and appear qualitatively 
correct, their practical value as predictive and quantita- 
tive tools is, at the present time, largely open to ques- 
tion. 

In summary, the difficulties hampering the practi- 
cal application of Navier-Stokes methods, including 
TRANSAM, are threefold: 

1) The lack of general accuracy and reliability of 
present turbulence models. 

2) The difficulty of generating the computational 
mesh or grid which engulfs the aircraft. 

3) The lack of adequate computer resources. 

FUTURE DEVELOPMENTS 
AND CONCLUDING REMARKS 

In its present state of capability, CFD provides 
large benefits in the design process, particularly when 
used judiciously and synergistically with the wind tun- 
nel. But the full realization of the great potential offered 
by the application of CFD to aircraft design poses for- 
midable, but not insurmountable, technical challenges. 
Undoubtedly, the NAS will be highly instrumental in 
advancing CFD technology the level needed to success- 
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fully meet them. Some of the major difficulties to be 
overcome, and probable future developments and 
trends are highlighted in the remainder of this paper. 

The most important obstacle to the ready applica- 
tion of advanced CFD methods to airplane design is the 
difficulty of generating adequate computational grids 
for complex three-dimensional configurations in a 
timely manner. Fortunately, significant progress is be- 
ing made in this area. Advanced graphics software and 
hardware developments, e.g., color graphics worksta- 
tions, are beginning to aid and speed up the grid gen- 
eration process. Finite- volume zonal methods, such as 
the ones discussed in this paper, facilitate the grid gen- 
eration task. In addition, alternate approaches to the 
treatment of complex geometries are being actively 
pursued, examples of which are the work discussed in 
References 11 and 12. 

Assuming continuing development, it is very likely 
that soon it will be possible to generate adequate com- 
putational grids about complete and arbitrarily complex 
aircraft configurations in a matter of hours or days rath- 
er than weeks ,or even months, which is the time scale 
of present grid generation capability. 

Turbulence modeling will continue to be the Achil- 
les' heel of Reynolds-averaged Navier-Stokes codes. 
Reference 13 provides a concise but comprehensive 
survey of the state-of-the-art in turbulence modeling. 
This survey makes obvious that many difficulties re- 
main. But in spite of the present short-comings of tur- 
bulence models, Navier-Stokes codes will play an in- 
creasingly useful role in many applications. 


Furthermore, turbulence modeling is one area where 
the NAS offers great potential for advancing the state- 
of-the-art. 

The availability of a wide variety of more capable 
scientific computers and supercomputers, coupled with 
advances in numerical solution algorithms, will greatly 
accelerate the application of Euler and Navier-Stokes 
methods to the solution of airplane design problems. 
Computers ranging in speed from 100 to 10,000 million 
instructions per second and in memory capacity from 
8 to 256 million words will become commercially avail- 
able at competitive prices in the near future. 

The analysis and optimization of multidisciplinary 
interactions will become commonplace thanks to the 
capabilities of CFD and the foreseen increases in com- 
puter power and cost-effectiveness. Coupling CFD 
with the computational methodologies of other disci- 
plines such as structures, propulsion, dynamics, and 
flight controls will allow the engineer to optimize not 
only the aerodynamics but also the structural compo- 
nents, propulsion and control systems, etc., in a fully 
interactive and integrated manner. Great performance 
and cost benefits are likely to result from the synergistic 
effects of this interactive multidisciplinary optimiza- 
tion. 

Finally, supercomputers and wind tunnels are 
complementary and not exclusionary tools. This has 
been repeated many times, but the more experience we 
gain with CFD applications, the more evident it be- 
comes. It is also true that to further advance CFD, more 
high-quality experimental data are required. 
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